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Presented are methods that avoid the need to employ an extrapolation technique in the region 
of the critical points for evaluation of the apparent emmissivity of diffuse cylindrical and conical 
cavities. The methods involve appropriate substitutions in the integrands of integral equations that 
are used in analytical solutions for determining the thermal radiation characteristics of diffuse and 
conical cavities. Equations for either isothermal or nonisothermal surface temperature conditions 
are provided in a direct form for computations. Numerical results are presented for a general linear 
temperature distribution along the length of a cylindrical cavity. The method is equally applicable 
for the solution of other problems in integral equations where discontinuities are encountered. 
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1. Introduction 

Analytical formulations have been derived [1, 2, 3] 1 
for determining the radiant interchange in finite length 
cylindrical cavities, and conical cavities whose bound- 
ing surfaces radiate in a gray diffuse manner. These 
solutions are in the form of integral equations where, 
unfortunately, most of the integrands exhibit discon- 
tinuities with finite limits at critical points, such as 
the corner of the cylindrical cavity and the apex of the 
conical cavity. Numerical solutions of these equa- 
tions in the region of the discontinuities have been 
previously evaluated by extrapolation techniques. 
Also, to a lesser degree of difficulty, some of the in- 
tegrands exhibit slope discontinuities, which for accu- 
rate numerical results may involve intricate numerical 
integration procedures. 

The purpose of this paper is to present methods for 
numerical evaluation of the apparent emissivity in 
these cavities that avoid use of an extrapolation tech- 
nique in the region of the critical points. For the main 
part, the methods involve appropriate substitutions 
which make the integrands go to zero at the critical 
points and at the points of slope discontinuity. Im- 
petus for this investigation has been the numerical 
evaluation of thermal radiation characteristics for 
parameters not included in the scope of references 1 
and 2. A particular need has been the determination 
of the thermal characteristics for low values of the 
surface emissivity of shallow cylindrical cavities in- 
vestigated by Kelly and Moore [4]. 



Also presented is a system of integral equations for 
determining the thermal characteristics of diffuse 
cylindrical and conical cavities whose surfaces are 
nonisothermal. A few selected numerical results are 
included for the cylindrical cavity with an arbitrary 
linear temperature variation over its length. 

2. Cylindrical Cavity 

Analytical formulations for the apparent emissivity 
of the surfaces of cylindrical enclosures have been 
derived by Sparrow, Albers, and Eckert [1]. Using 
the same nomenclature, their eqs (8) and (9) are pre- 
sented here in a different form. 

[L/d 
Ca(Xo) = € + ( 1 — €) 6 a (x)Ki (x , X)dx 

H — ^— {-,— Xo\ \ € a (r)K 2 (xo> r)rdr (1) 

1 — e f L ' d (L \ 

€ a (r) = e + — — I € a (x)K 2 (x, r) I — — x I dx 

, V( x , \x — x\ 2(x — x) 2 -\-3 
where Kfa> 9 x)=l — {{x(j _ xf+ 1}3/2 



(2) 



K 2 (x,r) = 



(z + w) 



{z 2 + {\-2w)z + w 2 } 



3/2 



1 Figures in brackets indicate the literature references at the end of this paper. 
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and e a is the apparent emissivity of a surface location; 
e a (xo) and e a (r) represent values on the curved surface 
and flat surface, respectively, of the cylindrical cavity. 
The physical system is described in [1]. L is the 
length of the cavity, d is the diameter of the cavity, 
x is a ratio of the distance measured along the cavity 
wall from the opening to the diameter, and r is the 
ratio of the radius of a position on the bottom of the 
cavity to the cavity radius, R. 

Initial attempts by the author at obtaining numerical 
solutions of (1) and (2) were frustrated by the erratic 
behavior of the kernels K 2 {xo, r) and K 2 (x, r) as r ap- 
proaches unity and x and Xo approach Ljd (the corner 
of the system). Visual inspection of these kernels 
and integrals involved indicates an indeterminate 
form at these limits. In correspondence with L. U. 
Albers (coauthor of [1]), it was disclosed that an 
extrapolation technique was employed in the region 
of the corner, and that the results were quite insen- 
sitive to wide variations in types of extrapolation. 
Attempts at extrapolation by the author proved to be 
inconclusive. Indications were that the second dif- 
ferences were quite large in this region, showing that 
linear extrapolation may be inaccurate. 

The suggested method for obtaining solutions is 
to make the integrand zero at the corner and make 
contributions to the results negligibly small in the 
region of the corner. This may be performed by 
making the substitutions 






(3) 

(4) 



in the integrals of (1) and (2). A , A u B , and B x are 
chosen such that 0(L/d) = </>(O) = and iftl) = 0(O) = O. 
Substitution of the first two terms of (3) and (4) in the 
appropriate integrals of (2) and (1), and by rdr = — 2dw, 

and \—j — x) dx = — dz/2 gives 



= 24 /2 f 

Jo 



^ (Bp + Biw)(zo + w)dw 
{w 2 -2z w J rzo + z 2 } 3 l 2 



I 2 = A + Aiw + 



[(L/d) 2 -w] [Ao + A^A^-AjL/d) 2 
[{(Lld) 2 -w} 2 + (L/d) 2 yi 2 



-h/^/2 In [2{[(L/d) 2 -w] 2 + (Lid) 2 yi 2 + 2(Lld) 2 +l-2w] 

(6) 

The kernel Ki(xo, x) has a discontinuous derivative 
at x = x . This did not seem to have an appreciable 
effect on the accuracy of numerical integration by 
Simpson's rule if the subdivision of L/dwas sufficiently 
small. The smaller the increment, the longer becomes 
the time needed to compute apparent emissivities, 
and it is therefore expedient to use as large an incre- 
ment as possible. This can be done by making the 
integrand equal to zero at x = Xq, and quite small in 
the vicinity of Xo, or by the substitution 

€ tt (x) = € a (Xo) + X(x) (7) 

in the first integral of eq (1). Letting 

r fxo cm I 

I:i = € a (xo)\ I K\(xo, x)dx+ \ K\(xiu x)dx\ 

and performing the indicated integration, we obtain 



h = €a(Xo) 



L 2xl+l 2 (l~*°) 2+1 

d 2(xf,+ l)" 2 „ f IT- 1»' : 



H-r 



With the substitution of (3), (4), and (7), eqs (1) and (2) 
become 



€a(x () ) =€ + (!—€) 



CM 

Jo 



\(x)Ki(xo, x)dx 



~r- f 7 i + (f -*<>) \l *MK»(*o, r)rdr 



4 
e„(r) = € + 



( -^ { /2 + / () ' " 9(x)K 2 (x, r)(~x 



(2a) 



r (L/d)2 



(Ap-i- Ajz)(z-\- w)dz 
{z 2 + (l-2w)z + w 2 }*l 2 



From a table of integrals 

(2z -h l)fl + 2z (zo -tMi 



/ I =2f 



Uo+1) 1 



-^/ 2 {2Bo+(l+22o)B,} 
-2i/ 2 Biln [2{(zo+D 1/2 - 



■»}]] 



(5) 



Numerical solutions were obtained from (la) and 
(2a) by the process of iteration, where all numerical 
integrations were performed by Simpson's rule and 
employed an increment equal to Ve4 of the difference 
of the limits of integration. This increment was 
proved adequate by comparison with results obtained 
using an increment of Vi28. For a fixed value of L/d, 
the kernel functions were computed only one time for 
all variations in the surface emissivity, €. Ki(xo, x) 
was computed in a one-dimensional array where proper 
indexing gave translation about a given value of Xo. 
K 2 (x, r) was computed in two-dimensional array. 

Initial values of e a (x) were determined from rather 
crude polynomial approximations. Placing the values 
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FIGURE 2. Apparent emissivity in a cylindrical cavity, L/d = l. 



FIGURE 1. Apparent emissivity in a cylindrical cavity, L/d — 2. 



for e a (x) in (2a), gave values for e a (r), which with the 
values of e a (x) placed in (la) gave new values for e a (x). 
This process was repeated using the new values in 
(2a). The number of iterations (rc) necessary for con- 
vergence was set by the criterion 



Table 1 



€ a {x=0)n — e a (x = 0)n- 
€ a (x = 0)„ 



; 0.0005 



which was usually satisfied for n = 4 or less. A typical 
time to compute 16 cases involving various values of 
L/d and e was about 2 min for an IBM 7094 digital 
computer. 

Figures 1, 2, 3, and 4 give values for apparent emis- 
sivity versus dimension ratios for L/<£=2, 1, 0.5, and 
0.25, respectively, and for various values of the surface 
emissivity, €. Where the same parameters occur, 
there is agreement with values given in table 1 of 
reference [1] to within 0.1 percent. Figure 5 is a plot 
of e a (r = 0) versus e for L/d = 2, 1, 0.5, and 0.25. Figure 
6 is a plot of the ratio of the overall emissivity of the 
cavity to the surface emissivity (QIlt€ttR 2 T a ) versus 
e for L/d= 2, 1, 0.5, and 0.25. The ratio was computed 
by numerical integration of eq (lib) of [1]. 









Values of €„(r = 0) and t„(r= 1) 






€ 


L/d-- 


= 2 


1.0 


0.5 


0.25 


0.9 


0.9936 


0.9944 


0.9785 0.9843 


0.9482 


0.9707 


0.9191 0.9600 


.8 


.9860 


.9877 


.9532 .9659 


.8924 


.9371 


.8361 .9149 


.7 


.9764 


.9795 


.9231 .9435 


.8314 


.8977 


.7506 .8633 


.6 


.9637 


.9686 


.8862 .9153 


.7640 


.8504 


.6619 .8030 


.5 


.9460 


.9534 


.8396 .8778 


.6879 


.7916 


.5694 .7316 


.4 


.9196 


.9302 


.7789 .8264 


.6005 


.7168 


.4719 .6452 


.3 


.8766 


.8913 


.6951 .7504 


.4976 


.6185 


.3684 .5388 


.25 


.8442 


.8612 


.6401 .6980 


.4383 


.5568 


.3138 .4758 


.2 


.7986 


.8179 


.5728 .6314 


.3724 


.4841 


.2570 .4046 


.15 


.7313 


.7526 


.4874 .5438 


.2982 


.3973 


.1977 .3237 


.1 


.6238 


.6457 


.3756 .4246 


.2136 


.2919 


.1354 .2311 


.05 


.4305 


.4488 


.2228 .2554 


.1156 


.1624 


.0697 .1242 



Values of QI(a€TrR 2 T 4 ) 



0.9 


1.0829 


1.0801 


1.0688 


1.0483 


.8 


1.1829 


1.1755 


1.1482 


1.1025 


.7 


1.3060 


1.2908 


1.2410 


1.1637 


.6 


1.4620 


1.4330 


1.3507 


1.2337 


.5 


1.6662 


1.6156 


1.4844 


1.3138 


.4 


1.9482 


1.8540 


1.6480 


1.4070 


.3 


2.3670 


2.1877 


1 .8543 


1:5163 


.25 


2.6609 


2.4072 


1.9788 


1 .5784 


.2 


3.0550 


2.6770 


2.1220 


1.6465 


.15 


3.6080 


3.0207 


2.2880 


1.7213 


.1 


4.4500 


3.4720 


2.4840 


i .8040 


.05 


5.8980 


4.0920 


2.7160 


1 .8960 



141 





1 1 
€=0.9 _____ 


1 1 




1 












\ 


~~~^ 


.9 






\ 
\ 


x^ 


.8 






\ 
\ 


\v 


.7 


0.4 




\ 




.6 

€ a (X) 
or 




€=Q3^^ > 
0.25^^^^^ 


\ 


!\ v ^ 


€ a (r)-5 




0.2^-""'" 


\ 


\j 


4 


- ^^ 


0.15^-^ 


\ 


V VI 


.3 


.1 




\ 


v x ^- 


.2 


- 


X>5___- — "~"~"~ 


\ 


^J 


.1 


1 1 


L/d =0.5 

1 1 




i 





1 1 1 1 


- 


1 




0^9 — 




_____ — 


\ 




.9 


- ^ — ■— " 


\ 




_____ ^ 


\ 






CX7____ — — - ~~ 


\ 




.8 


06__^— ^^~^ 


\ 




.7 




\ 




.6 






\ \. 


6 a (X) 




\ 




or 


03^-^^ 






6 a (r) -5 


025__^^^ 


\ 

\ 


\\^ 


.4 


(X2_^^ 


X 


\\^ 


.3 


- J5__^^^^^^ 




Vxj 


.2 


^ — ^^ 






.1 


L/d=.25 

1 1 1 1 








i 



4 .6 .8 

X/L 



IJO 1.0 0.5 

r/R 



.6 
X/L 



.8 1.0 1.0 



0.5 
r/R 



FIGURE 3. Apparent emissivity in a cylindrical cavity, L/d = 0.5. FIGURE 4. Apparent emissivity in a cylindrical cavity, L/d = 0.25. 
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FIGURE 5. Apparent emissivity versus surface emissivity at the center 
of the bottom of cylindrical cavity. 
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FIGURE 6. Overall heat transfer rate from cylindrical cavities versus 
surface emissivity. 



3. Conical Cavity 

Analytical solutions for the apparent emissivity 
of the surfaces of diffuse conical cavities were de- 
rived by Sparrow and Jonsson [2], Here, the problem 
is to evaluate an integral of the form 



(0/2) 



2y sin (6/2) J„ 



e (l (t)Kz(y, t)dt 



(9) 



where 



K,(y,t)=l-\t-y\ 



[(*-y) 2 + 6fysin 2 (fl/2)] 
[(t-y) 2 + 4>tysm 2 (6l2)fl 2 



and y is the ratio of a surface position measured from 
the apex to the total length of a cone. The kernel 
function K^y, t) has a discontinuous derivative at 
y= t and a problem arises for the case y=0, which by 
visual inspection gives an indeterminant form. This 
suggests a substitution of the form 



e„(t) = €„(y) + y{t). 



(10) 



Letting 



\" K-ly, t)dt+ \ K 

J() J if 



s(y, t)dt 



2y sin (012) 

and performing the indicated integration, we obtain 
2 sin (0/2) -1 , cos 2 (0/2) 



h(y)-- 



2 sin (612) 2y sin (6/2) 

(l-y) 2 + y(l-fy) tan 2 (fl/2) 



1 



{(l-y) 2 + 4ysin 2 (<9/ 2 )} 1/2 



(11) 



For the case y = 0, the second term is an indeter- 
minant, but can be evaluated by L'Hospital's rule. 
This yields 



/.,(())= i-sin 3 (8)2). 



The integral (9) then becomes 



e„(y)h(y) - 



cos 2 {012) 



2y sin (0/2) 



|/>*3(y, 



(12) 



t)dt (13) 



where the second term still contains the indetermi- 
nant form for y = 0. By letting t > 0, and applying 
L'Hospital's rule, we find that the integrand goes to 
zero. Also, y(0) = 0, so that the second term goes to 
zero for y=0, and the limit becomes 



Lim cos 2 (6/2) 
y-»0 2ysin (0/2) 



f 6 " 



(t)K,(y, t)dt ■■ 



:e„(O)[l-sin 3 (0/2)]. 
(14) 



This is the same result derived by C. H. Page in an 
internal NBS Report (1952), where it was intuitively 
assumed that the apparent emissivity at the apex of 
the cone would be independent of the length of the 
cone, or that the cone appears infinite. The upper 
limit of (9) was changed to infinity; this yields e„(t) 
constant, hence equal to e (l (0). Another way to arrive 
at this limit for the case y = is to assume y(t) = A\t 
-\-A 2 t 2 + . . . in eq (13) and perform the indicated 
integration and limiting process which will give the 
same result as eq (14). 

By making the substitution y = x/L in eq (5) of [2], 
the equation for the apparent emissivity of a cavity 
maintained at uniform temperature is 



e„(y) = € + (l-€) \e (l (y)I 4 (y) 
cos 2 (0)2) 



+ 



2y sin (0j2) 



12) Jo 



y(t)K,(y, t)dt\ (15) 



which has a unique solution at y=0, 
e«(0) = 



e + (l-e) sin 3 (0/2) 



(16) 



This is the expression for conical cavities shown by 
F. J. Kelly, in a paper soon to be published, as derivable 
from the expression for cavities used by Gouffe. It 
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is also the expression derived by C. H. Page by another for various apex angles is given in figure 7, where the 
method (see discussion following eq (14)). Com- average deviation over the length is — 2.3 percent for 
parison with graphical quantities given in [2] for e = 0.3 and = 45°. 
e[l — e„(0)]/(l — e) show almost identical agreement. 
As in [2], numerical solutions were obtained by the 
process of iteration, where initial values of e„(y) were 4 Cavities With Nonisothermal Surfaces 
computed irom 

Equations 1, 2, and 15, are concerned with cavities 
(17) for which the surfaces were assumed isothermal. 
In the succeeding paragraphs, there are presented 
systems for determining the apparent emissivities of 
cylindrical and conical cavities with arbitrary varia- 
tions of surface temperature in respect to dimensions 



e " (v) ' l-(l-€)/ 4 (y)* 
Because graphical presentations have appeared in 



[2], numerical values are not shown in this paper. 

:_ Eq . Uati ^n ( rl 7 ]_ Cai ]i be C ™ sid Jl d _ * J°° d 5i!?™^; oTTength 1^^!^^ a7sumlng the7u7face^m7^ 

is constant and not a function of temperature. 



tion to eq (15) for all cavities with surface emissivities 
equal to or greater than 0.7 and for all cavities with 
apex angles greater than 120° and emissivity greater 
than 0.3. For emissivities and apex angles below 
the values cited above, the deviation between the two 
equations is not acceptable. An example of the per- 
centage deviation 



y(y) z 



100 {ea{y)-ea{y),~\ 
£ a {y) 



(18) 



4.1 Cylindrical Cavity 

For the cylindrical cavity, the integral equations 
take the form 



B(xo) 



e a (xo) = -pfi = *fM + (1 - €)Ci [e«M, e a (r)] 



(19) 
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FIGURE 7. Percentage deviation of eq (17) from eq (75) versus position 
from apex of a conical cavity. 

Unless otherwise noted, surface emissivity. e = 0.3. 



e«(r) = ^j=eg(r) + (l-e)G 2 [e a (*)] 



(20) 



where 



T*(x ) = nf(xo):f(xo) = Co + C iy +. . . + C n y m (21) 

r(r) = rte(r);^)-l + D^+. . .+D k r", (22) 

y = Xol(Lld), To is an arbitrary temperature which is 
taken to be the temperature at r = 0, and G\ and G 2 
are the coefficients of (l-e) in (la) and (2a), respec- 
tively. Two sets of integral equations are defined by 

E n (xo) = *r + (1 - e)Gi [£„(*), E n (r)] (23) 

E n (r) = e + (l-e)G 2 [E n (x)] (24) 

n = Q, 1, 2, . . . m 

Mj(xo) = e + (1 - e)Gi[M/(*), M{r)] (25) 

Mj(r) = eyV + (1 - e)G 2 [Mj{x)\ (26) 

7=1,2, . . . k 



and 



The apparent emissivity is then found from the 
following: 

e a (xo) = CyEtJixo) 4- &Ei(xq) + . . . + C m E m (x ) 

+ DiM l (x { >) + . . .+D k M k (x {) ) (27) 

ta(r) = C; } E {) (r)+C l E l (r) + . . . + C m Ejr) 

+ D t Mi(r)+. . . + D k M k (r) (28) 
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E n (X ) 




i.o 



X/L or r/R 

FIGURE 8. Component apparent emissivity values versus position in 
cylindrical cavity (L/d = 2.0) for a surface emissivity = 0.9. 

which satisfies (19) and (20) for C' Q = C Q -Dx-D 2 - . . . 
Dk, which gives continuity to the temperature at 
x { ) = L/d and r=l. Figures 8 and 9 show solutions 
of (23) and (24) for € = 0.9 and 0.5, respectively, with 
L/rf=2.0, and n = 0, 1, 2, 3, and 4. 

One simplifying procedure is to assume that the 
bottom of the cavity is isothermal at temperature T () . 
This eliminates the numerical solutions of (25) and 
(26); and (27) and (28) become 
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FIGURE 9. Component apparent emissivity values versus position in 
cylindrical cavity (LJd = 2.0)for a surface emissivity =0.5. 

The ratio of the cavity radiant heat exchange Q with 
the environment through the opening to that of a black 
disk of equal area at temperature T {) is 



o-7>K 2 (1-6) 



I n 



(xo)}dx 



+ 2 I* {l-e a (r)}rdr 



£C„<? H (31) 



e«(* ( >) = ^ CuEdxt!) 



e„(r)=2c„£„(r). 



(29) 



where 



(30) <?„ = 



1-e 



4^ f ' {%8-£„(*„)}Jx+2 | °{l-.E„(r)}rdr 

a Jo Jo 
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TABLE 2. Component thermal characteristics for a cylindrical 
cavity, h/d — 0.5, with a nonisothermal wall 





6 = 


0.5 


0.6 


0.7 


0.8 


0.9 




n = 


0.68786 


0.76398 


0.83139 


0.89236 


0.94823 




1 


.59722 


.68197 


.76344 


.84298 


.92158 


E„(r = 0) 


2 


.56747 


.65522 


.74133 


.82695 


.91295 




3 


.55334 


.64260 


.73095 


.81946 


.90893 




4 


.54539 


.63554 


.72516 


.81530 


.90670 




n = 


.79159 


.85038 


.89769 


.93707 


.97066 




1 


.74185 


.80906 


.86620 


.91592 


.96006 


E„\r=\) 


2 


.72380 


.79416- 


..85485 


.90831 


.95627 




3 


.71433 


.78637 


.84894 


.90436 


.95421 




4 


.70850 


.78158 


.84530 


.90193 


.95310 




n = 


.68571 


.76237 


.83139 


.89188 


.94810 




1 


.11781 


.10232 


.08175 


.05731 


.02986 


E„U,. = 0) 


2 


.09700 


.08423 


.06729 


.04718 


.02459 




3 


.08713 


.07574 


.06057 


.04252 


.02219 




4 


.08143 


.07087 


.05674 


.03988 


.02084 




n = () 


.7422 


.8104 


.8687 


.9186 


.9619 




1 


.4802 


.5223 


.5576 


.5876 


.6136 


<?. 


2 


.4006 


.4362 


.4662 


.4921 


.5149 




3 


.3630 


.3959 


.4240 


.4484 


.4703 




4 


.3412 


.3728 


.3999 


.4238 


.4453 



Tables 2, 3, and 4 give values of the component 
thermal radiation characteristics E„(r= 0), E,,(r= 1), 
E H (xo=0) and Q„ for L/d=2. 1. and 0.5, 6 = 0.5, 0.6, 
0.7, 0.8, and 0.9, and n = Q. 1. 2, 3. and 4. For the 
sake of a numerical example, assume a cylindrical 
cavity, € = 0.7, where the temperature decreases 
linearly from a temperature T {) at its base to a value 
0.97 7 ,, at the open end, or, from (21), /(*) = 0.6561 + 
0.2916y+ 0.0486/ + 0.0036/ + O.OOOly 4 . For L/d=2, 
from eq (31) and table 4, we find 

Q 



rrTlrrR 2 



= 0.6561 x 0.9142 + 0.2916 x 0.2835 



0.048 X 0.1628 + 0.0036 X 0.1210 



+ 0.0001X0.1012 = 0.6908 

which is compared to a value of 0.9142 for the iso- 
thermal cavity at 7V 



Table 3. 



Component thermal characteristics for a cylindrical 
cavity, L/d = 7, with a nonisothermal wall 





e = 


0.5 


0.6 


0.7 


0.8 


0.9 




ra = 


0.83957 


0.88625 


0.92309 


0.95324 


0.97847 




1 


.70223 


.76827 


.82950 


.88782 


.94438 


£„(/• = 0) 


2 


.64592 


.71951 


.79050 


.86032 


.92993 




3 


.61520 


.69291 


.76914 


.84521 


.92196 




4 


.59658 


.67656 


.75593 


.83583 


.91700 




n = 


.87784 


.91531 


.94353 


.96591 


.98431 




1 


.79987 


.85305 


.89741 


.93563 


.96942 


E„['r=l) 


2 


.76831 


.82782 


.87869 


.92333 


.96338 




3 


.75077 


.81384 


.86831 


.91651 


.96003 




4 


.73989 


.80500 


.86170 


.91216 


.95789 




ra = 


.70148 


.77174 


.83533 


.89394 


.94858 




I 


.09719 


.08009 


.06114 


.04117 


.02068 


£„Un = 0) 


2 


.06924 


.05627 


.04237 


.02816 


.01396 




3 


.05704 


.04613 


.03455 


.02284 


.01127 




4 


.05059 


.04075 


.03043 


.02007 


.00988 




n = 


.8078 


.8598 


.9036 


.9404 


.9721 




1 


.4088 


.4235 


.4334 


.4401 


.4447 


Q„ 


2 


.3018 


.3100 


.3149 


.3176 


.3189 




3 


.2552 


.2614 


.2651 


.2671 


.2682 




4 


.2281 


.2343 


.2382 


.2405 


.2419 



4.2. Conical Cavity 

For the conical cavity, the integral equation takes 
the form 

ea(y) = ^- = ef(y) +(1 - e)G,[e (1 (t)} (32) 

where 

Hy) = Ttf(y) = 1 + Cry + C 2 f + . . . + C m y m (33) 

and To is the temperature at the apex of the cone, 
Gslead)] is the coefficient of (1 — e) in (15) and B(y) is the 
radiant emission. 

A set of m integral equations is defined by 



£„(y) = e/>+(l-e)G 3 [e a W] 
/i = 0, 1, 2, . . ., m. 



(34) 



Table 4. Component thermal characteristics for a cylindrical 
cavity, L/d = 2, with a nonisothermal wall 







e = 


0.5 


0.6 


0.7 


0.8 


0.9 




n = 


0.94601 


0.96371 


0.97639 


0.98599 


0.99364 




1 


.81188 


.85575 


.89514 


.93171 


.96644 


Ejr = 0) 


2 


.74333 


.79923 


.85163 


.90204 


.95130 




3 


.70026 


.76313 


.82354 


.88269 


.94132 




4 


.67076 


.73806 


.80386 


.86902 


.93423 




n = 


.95338 


.96864 


.97949 


.98773 


.99436 




1 


.86738 


.90331 


.93281 


.95796 


.98006 


£>=1) 


2 


.82589 


.87139 


.90974 


.94312 


.97288 




3 


.80054 


.85166 


.89543 


.93388 


.96840 




4 


.78341 


.83819 


.88559 


.92750 


, .96530 




n = 


.70647 


.77428 


.83653 


.89437 


.94868 




1 


.06475 


.05108 


.03755 


.02455 


.01204 


£«<*„ = 0) 


2 


.03595 


.02713 


.01920 


.01209 


.00573 




3 


.02561 


.01890 


.01309 


.00809 


.00376 




4 


.02066 


.01500 


.01034 


.00633 


.00292 




/( = () 


.8331 


.8772 


.9142 


.9463 


.9746 




1 


.2884 


.2863 


.2835 


.2801 


.2766 


Q« 


2 


.1737 


.1679 


.1628 


.1580 


.1534 




3 


.1309 


.1258 


.1210 


.1167 


.1130 




4 


.1087 


.1051 


.1012 


.0977 


.0948 



Then the apparent emissivity is 

€a(y) = Eo(y) + C.E^y) + . . . + C m E m {y). (35) 

Evaluation of (34) for substitution in (35) gives an exact 
solution to (32). An approximation from eq (17) useful 
within the limits cited in section 3 is 



e«(y) = 



ef(y) 



l-(l-e)/ 4 (y) 



(36) 



5. Discussion 

Sections 2 and 3 present an analytical treatment by 
appropriate substitution in integrals which exhibit 
slope discontinuities or apparent discontinuities or 
both at critical points, such as the corner of a cylin- 
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drical cavity and the apex of a conical cavity. The 
substitutions (3), (4), (7), and (10) make the appropriate 
integrands go to zero when evaluated at the critical 
points and points of slope discontinuity, whereby the 
awkward behavior of the integrand is practically elimi- 
nated. With the transformed problem, as given by 
(la), (2a), and (15), accurate numerical results are 
readily obtained by use of a sufficiently small integra- 
tion step and the fulfillment of the convergence cri- 
terion for the successive iterations. 

It may be remarked that the very good agreement 
found between the numerical results arrived at by 
extrapolation in references 1, 2, and 3, and the rigorous 
analytical formulation given here shows that any error 
introduced by linear extrapolation was negligible. 
However, to avoid the uncertainties involved in ex- 
trapolation, this paper does provide the equations 
(isothermal and nonisothermal) in a direct form for 
computation. The suggested method is equally ap- 
plicable in other cases of integral equations in which 
similar discontinuities are encountered. 

The nonisothermal cylindrical cavity has been 
treated by Sparrow [5]. For the same values of the 
parameters, there is generally good agreement between 
the numerical results presented in this paper and 
reference 5. Two assumptions were made in the 
treatment [5]; namely, (1) the temperature over the 
base of the cavity was isothermal, and (2) the apparent 
emissivity over the base was constant. For relatively 



deep cavities, these assumptions can be considered 
reasonable, but for shallow cavities there is some ques- 
tion as to the validity of these assumptions. Figures 
2, 3, and 4, for example, show a considerable variation 
of apparent emissivity e a (r) over the base of the cavity. 
Although numerical results are not presented in this 
paper for temperature variations over the base of the 
cavity, this can be accomplished by evaluation of eqs 
(25) and (26). Also, numerical results are presented 
for general linear temperature distributions, but may 
be evaluated for all temperature distributions that may 
be represented in a polynomial or transcendental form. 
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